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Abstract 

In this paper, via the generalized Darboux transformation, rational soliton solutions are derived for 
the parity-time-symmetric nonlocal nonlinear Schrodinger (NLS) model with the defocusing-type non¬ 
linearity. We find that the first-order solution can exhibit the elastic interactions of rational antidark- 
antidark, dark-antidark, and antidark-dark soliton pairs on a continuous wave background, but there is 
no phase shift for the interacting solitons. Also, we discuss the degenerate case in which only one rational 
dark or antidark soliton survives. Moreover, we reveal that the second-order rational solution displays 
the interactions between two solitons with combined-peak-valley structures in the near-field regions, but 
each interacting soliton vanishes or evolves into a rational dark or antidark soliton as \z\ oc. In 
addition, we numerically examine the stability of the first- and second-order rational soliton solutions. 
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1 Introduction 


In classical quantum mechanics, a basic assumption is that the Hamiltonian operator is Hermitian to ensure 
that every physical observable is associated with a real spectrum [1]. In 1998, Bender and Boettcher proved 
that a non-Hermitian Hamiltonian also has real and positive eigenvaluesprovided that it has the combined 


parity and time reversal symmetry (usually called the VT symmetry) [h]. Such pioneering work has led 


to the complex extension of quantum mechanics [3]. In general, a necessary condition for a Hamiltonian 


H = p 2 /2 + V(x ) to be PT-symmetric is that the complex potential satisfies V(x) = V*(—x), where 

Also, the notion of VT symmetry has been applied to other 


p denotes the momentum operator 
areas of theoretical physics, including Lie algebra (5^, complex crystals (fj, quantum chromodynamics 




equation in quantum mechanics 
ing and testing the VT -related concepts 


n 


&, 

Bose-Einstein condensates (8[], classical mechanics (9t], and so forth. 

Owing to the similarity between the paraxial equation of diffraction in optics and the linear Schrodinger 

ildlllh. it is regarded that optics can provide a fertile ground for realiz- 
An optical potential respecting the VT symmetry can 
be realized in the complex refractive index distribution no + ur(x) + in/(x), where the real index pro¬ 
file ur(x) should be an even function while the gain or loss component n/(x) must be odd, and no is a 


constant background index 


HQ. it 


has been shown that VT -symmetric optical structures can exhibit 


unique characteristics such as double refraction, power oscillations, spontaneous VT symmetry breaking, 


nonreciprocal diffraction patterns, and unidirectional invisibility [11-15]. In experiments, VT symmetry 
breaking within the realm of optics has been observed Il3l. which has stimulated the development of 
VT optical materials and optical elements [16, 17]. 

In nonlinear optics, the VT symmetry has received considerable attention in the last few years 111814301. 
Musslimani et al. first suggested the existence of optical solitons with the presence of the Scarff II poten¬ 
tial and periodic PT-symmetric potential [18]. Later on, a lot of work was devoted to the existence and 
stability of nonlinear modes in different VT - symmetric nonlinear systems,_such as the fundamental and 
higher-order solitons in the Gauss [19], harmonic [20], and Rosen-Morse 12ill VT- symmetric potentials, 
gap solitons |22] and defect solitons [23] in periodic PT-symmetric potentials, localized modes supported 
by VT -symmetric nonlinear lattices 124], lattice solitons in VT - symmetric mixed linear-nonlinear optical 
lattices [25], and vector solitons [26], breathers [27] and rogue waves [28] in PT-symmetric coupled waveg¬ 
uides. Meanwhile, researchers have studied the effects of nonlinearity on the VT symmetry breaking [29] 
and dynamical characteristics of a beam in PT-symmetric optical nonlinear systems 13011 . 

Recently, Ablowitz and Musslimani proposed the following VT -symmetric nonlocal nonlinear Schrodinger 
(NLS) equation [31]: 


i u z (x, z) = u xx (x , z) + 2 eu(x, z)u*(—x, z)u(x , z) (e = ±1), (1) 

which is obtained from the standard NLS equation by replacing \u\ 2 u with u(x, z)u*(—x, z)u(x , z), where 
u(x,z) denotes the electric field envelope, z is the spatial coordinate along the propagation axis, x is 
the transverse coordinate, e = ±1 denotes the focusing (+) and defocusing (—) nonlinearity, and the 
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star signifies the complex conjugate. The VT symmetry of Eq. A]) means that the self-induced potential 
V(x,z) = 2 eu(x, z)u*(—x, z) exactly satisfies the relation V(x,z) = V*(—x,z), while the nonlocality 
says that the value of the potential V(x, z) at x requires the information on u(x, z) at x as well as at — x Q. 
Equation ([!]) is integrable in the sense of admitting the Lax pair and an infinite number of conserved quan¬ 
tities. Therefore, its initial-value problem can be solved by the inverse scattering transform (1ST) 113ill . In 
addition, the integrability of the discrete version of Eq. ([!]) has been established [32, 33], and some other 
PT-symmetric nonlocal integrable models have also been proposed [34]. 

In contrast with the standard NLS equation, the VT -symmetric nonlocal NLS model has many different 
properties. For the focusing case, Eq. (JT|) possesses both static bright and dark solitons 132 l I35H. but its 
moving soliton obtained by the 1ST method contains a singularity Dill . In the defocusing case, via the 
elementary Darboux transformation (DT), we have revealed the nonsingular exponential soliton solutions 
on a continuous wave (cw) background [36]. Such exponential solitons in general appear in dark-dark, 
antidark-antidark, antidark-dark or dark-antidark pairs, and can exhibit the usual elastic interactions. It 
should be pointed out that the antidark soliton is another type of soliton existing in the normal dispersion 
regime on a nonzero cw background [37]. However, the exponential solitons in Eq. (OD) will become unstable 
if the solution has a small shift from the center of the VT symmetry, which has been confirmed by numerical 
simulation 132, 36]. 

We note that the kernel for constructing the elementary DT is made up of linearly independent solutions 
of the Lax pair associated with different spectral parameters [38]. Hence, the elementary DT cannot deal 
with the degenerate cases when the spectral parameter in the Lax pair reduces to some fixed value, so that 
the explicit solutions are derived only in the exponential form (similar to the soliton and breather solutions). 
In fact, via the generalized DT proposed by Matveev [39U, one can obtain the rational solutions in such 
degenerate cases. In recent years, the generalized DT has been widely used to construct the rational rogue 
wave solutions (which are algebraically localized in any direction of the temporal-spatial plane) of the NLS- 
type models ll4(i klh. In this paper, we will construct the generalized DT of Eq. ([I]) based on the work in 
Ref. [36], and further reveal the rational soliton phenomena on the cw background. In sharp contrast with 
the rogue wave solutions, the rational solutions obtained in this work are localized along the straight lines in 
the xz- plane, and can display the profiles of the common exponential dark and antidark solitons. Therefore, 
we call these two types of localized wave structures as the rational dark (RD) and rational antidark (RAD) 
solitons, respectively. Different from the exponential soliton solutions obtained in Ref. [36], the first-order 
rational soliton solution exhibits only the elastic interactions for the RAD-RAD, RD-RAD, or RAD-RD 
soliton pairs, but the interacting solitons do not experience the phase shift. We also find that the second- 
order rational solution displays the interactions between two solitons with combined-peak-valley structures 
in the near-field regions, but each interacting soliton will eventually vanish or evolve into a RD or RAD 
soliton as \z\ —»• oo. Our numerical experiments show that the rational soliton solutions have good stability 
against small initial perturbations, but their stability will be destroyed if the self-induced potential loses the 
VT symmetry with respect to any point of x. 

The structure of this paper is as follows: In Sect. [2l we will construct the generalized DT based on the 
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elementary one in Ref. L3£j]. In Sect.O we will derive the first- and second-order rational soliton solutions 
of Eq. (0D) with e = — 1, discuss the soliton interaction properties via asymptotic analysis, and examine the 
stability of the rational soliton solutions by performing numerical experiments. In Sect. ® we will conclude 
this paper. 


2 Generalized Darboux Transformation 


The Lax pair of Eq. (0D) can be written in the form 113 ill 

t ttt ( ^ u(x,z)\ 

VE r = [Ap = ' 1 




-£U*( — x,z) —A 

^ = v^> = ( -2iA2 ~ ieu ( x ’ z ) u *(~ x ’ z ) -2iA u(x,z)-iu x (x,z)\y 

\.2is\u*(—x,z)—isu*(—x,z) 2iA 2 + i eu(x, z)u* (—x, z) I ’ 


(2a) 

(2b) 


where V I' = ( f,g) T (the superscript T represents the vector transpose) is the vector eigenfunction, A is the 
spectral parameter, and Eq. (OQ) can be recovered from the compatibility condition U z — V X + UV — VU = 0. 

On the basis of work described in Ref. IBal . the Ath iterated elementary DT for Eq. Q} can be constituted 
by the eigenfunction transformation 


V[N]=T[N]V, T[N] = 


/ N N \ 

X N - E On(*,*)A "- 1 - E b n {x,z){-\Y~ l \ 




n =1 n=l 

N N 

- E C n (x,z)\ n ~ 1 A^ — E d n (x,z)(- A ) n_1 

n=l n=l > 


(3) 


and the potential transformation 

u[N](x,z) = u(x,z) + 2 (— u*[N](—x,z) = u*(—x,z) + 2 scj\r, 


(4) 


where N represents the iterated time. The new eigenfunction [N] is required to satisfy the Lax pair in 
Eqs. (l2al) and (l2bl) with u[N](x,z) and u*[N](—x,z) instead of u(x,z) and u*(—x,z), respectively. The 
functions a n (x, z),b n (x, z),c n (x, z), and d n (x, z) (1 < n < N) can be determined from 


T[N] | A=Afe = o, T[N] |a=a: = o (1 <k<N), 


(5) 


where 'Iq. = [fk{x, z), giJx, z)] T and 'Iq. = [gj,(—x,z),e /£ T are the solutions of the Lax pair in 


Eqs. (l2al) and (I2b1) with A = A/,, and A = A)’, respectively. In particular, via Cramer’s rule, the functions 
bpi(x, z) and cn(x, z) can be obtained in the determinant form 


h/v = (—1) 


N—l TN+1,N—1 

t~n,n 


Cn = 


tn- 1,N+1 
t~n,n 


with 


tn,n 


Fn xN Gn xJV 

zGnxN F Nx n 


( 6 ) 


(7) 











where the block matrices F NxN = [AJ* 1 fk(x, 2 )] 1<t G Kx n = [(-A0 m 'gt(x, 2 )] Kk m<K , 

Note that the elementary DT does not apply to the degenerate case when A^ 1+ i,..., Xk 2 -i X kl , 
Afc 2 +i5 • • • 5 A^ 3 _i — > X k2 , * * •, Afc n _j_i,..., Xn —>■ A& n , where 1 = Aq < &2 < • • • < k n < N, 1 < n < N, 


and X ki A&. (1 < i < j < ri). For such a degenerate case, the functions a n (x, z ), b n (x, z ), c n (x, z), and 
d n (x,z ) (1 < n < iV) in the Darboux matrix T[iV] cannot be uniquely determined because the coefficient 
matrix in Eq. ([5]) is singular. To overcome this problem, we define mi = ki+\ — k{ — 1 for 1 < z < n — 1 
and m n = N - k n , and assume that f k .+ h (x,z) = f ki {x,z, X ki+h ) and g ki+h (x,z) = g ki (x,z,X ki + h ) 
if mi > 0, where X ki+h = X ki + e*, 1 < h < m i9 e* are small parameters, and [f k .(x, z),g ki (x, z)] T 
corresponds to the solution of the Lax pair in Eqs. (l2al) and (l2bl) with X = X ki (1 < i < n). 

Via the idea of Matveev’s generalized DT |39j,[40], we expand the elements in the block matrices Fjy X N, 
GnxN, FnxN, and Gn x n in the Taylor series form 


A % + 1 h fk i +h(x,z,\ ki +h) = J2fk? > 

i=o 


(-A fci+/l ) m 1 <?fc i+ h(a;, 2 ,A fc . +fc ) = 5^^g* 1,J \x,z)ej, 

j =0 


(8a) 

(8b) 


with 


Am-lJ) 

J ki 


(x,z) 


(m— 1 , 7 ')/ \ 

ffL 


1 9 j [A^ + j t / fc .(x,z,A fc . +/t )] 

•? ! ei=0 ’ 

1 d j [(~A fc i+fe ) m ~ 1 gfc,(a;,^,A fc . +/t )] 


where l<i<n, l<m<./V, l</r< m,, and particulary fj^’°\x,z) = f^x^z), g^’°\x,z) = 
9ki(x,z). 

Taking the limit e* —>> 0, one can find that the coefficient matrix in Eq. (0) is no longer singular, which 
means that all the undetermined functions in T[N] can be uniquely solved. Thus, the new potential trans¬ 
formations are given as 


with 


where 


u[N](x, z) = u(x, z) + 2 N+ ^ N 1 ? u*[N](—x, z) = u*(—x, z) + 2s iV / 


iV-l,JV+l 


1 N,N 


1 N,N 


r N,N ~ 


r NxN ^NxN 
£ G'nxN F'nxN 



( F[^ 


( G \ N 


( Fp 



F’nxN — 


> g' N xN — 


1 F'nxN — 


5 G'nxN — 



UJ 


\G' n ) 


vJ 


\Gn/ 


(9) 


( 10 ) 


( 11 ) 
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and Fj, G[, F[, and G\ (1 < i < n) are the block matrices F[ = (x, z) 


0<p<iV 

0<q<m, 


^G' = 


F' = 


(- 1 ) P /fcf P ’ 5) (-*» 2: )J o< p <jv-i’ and ^ = [(- 1 ) P ^fe i U, ’ ,; (-*.«) 

0<g<m^ 


*(p,q)/ 


0<p<N-l ' 
0<g<m^ 



(x,z) 


0<p<iV — 1 ’ 
0<g<m^ 


Therefore, the eigenfunction transformation in Eq. ([3]) and potential transformation in Eq. ([9]) constitute 
the generalized DT for Eq. (0Q) when some of the spectral parameters {A ^}^ 1 coincide with each other. It 
is obvious that the elementary transformation in Eq. (0) corresponds to the particular case of the generalized 
one in Eq. © when n = N. To avoid the triviality of the DT, we require that A& cannot be a real number. 
In the next section, we will use the generalized DT to construct the rational soliton solutions of Eq. ([1} with 
£ = — 1 on a cw background. 


3 Rational Solitons on a cw Background 


It is easy to see that Eq. A]) with e 


-1 has the plane wave solution [36] 


u = pe Kx+[nz ^ (Q = 2 p 2 - K 2 ), 


( 12 ) 


where K, p, and (j) are three real parameters. Obviously, such a plane wave solution is unstable when 
K ^ 0. Thus, we take K = 0 in Eq. (fl2l) and substitute it into the Lax pair in Eq. d2al) and (l2bl) with A = A& 
(1 < k < N), which gives the solution for /& and as follows: 

e 2 (a k e^ kXk + /3 k e~ flkXk ) 

c _ 2i P z+i<j> ^ k -\ k )ak c Hk/Xk _ (Mfc+Afc)fifc c -/ij,Xk^ 

with p k = X\ + p 2 and Xk = x — 2i X k z, where a k and /?*. (1 < k < N) are free complex parameters. 
It should be noted that the solution in Eq. (ITT!) will reduce to a rational one if A& = icrp. Hence, in order to 
derive the rational solutions, we take Ai = lap and let A 2 ,..., Ajv all degenerate to Ai, which corresponds 
to n = 1, k\ = 1, and m\ = N — 1. Following the expansions in Eqs. (l8al) and (l8bl) . we can obtain the 
formulas for z) and z) (1 < m < N + 1; j = 0,1,2,...). For simplicity, in the 

following calculations, we set e* = S 2 , a k = ^ =1 1} , [3 k = —e~ flk ^ =1 1} for 1 < fc < iV, 

where Sj are arbitrary complex numbers. 




3.1. First-order rational soliton solution 

By truncating the expansions in Eqs. (l 8 al) and (l 8 bl) at j = 1, we obtain the first-order rational soliton 
solution as follows: 


(2p£ + iT — icr) (2pry — iT* + icr) 

+ 1(|K|2 + 1) 


(14) 


where £ = x + 2 crpz, 77 = x — 2 apz, K = 2psi + icr, and cr = ±1. It can be proved that the solution in 
Eq. (fl4l) has no singularity if and only if the parameter s\ satisfies the condition 

Im(si) ^ (15) 

2 P 
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Under this condition, we perform an asymptotic analysis of the solution in Eq. (fl4l) so as to clarify the 
dynamical behavior underlying the solution. 




(a) (b) (c) 


Figure 1: Three types of elastic two-soliton interactions via the solution in Eq. (fl4l) : (a) RAD-RAD interaction with 
p = 0.5, (j) = 0, a = 1, si = 1 — 0.2 i. (b) RAD-RD interaction with p = 1, </> = 0, a = 1, si = 1 — 2 i. (c) RD-RAD 
interaction with p = 1, 0 = 0, a = 1, si = 1 + i. 


First, we obtain the asymptotic expression of the solution in Eq. (fl4l) along the line x + 2 apz ^ 0 as 
z\ -A oo as follows: 


Ui 


M 2 


«i 




pe 


1 - 


lA 2 


1 - 


Ap(x + 2 crpz) + 4p si 
2p(x + 2 crpz) + K J 5 

8crplm(si) 


|2p(x + 2apz) + K | 2 


(16a) 

(16b) 


For the cases crlm(si) > 0 and crlm(si) < 0, the intensity \u\\ 2 can respectively exhibit the RD soliton 
beneath the cw background u = and the RAD soliton on top of the same background, and the 

valley and peak are both localized along the line x+2crpz+~Re(si) = 0. Along the line x J r2apz+Re(si) = 
0, | u\ | 2 reaches the minimum for crlm(si) > 0 and the maximum for crlm(si) < 0. Thus, the height of the 
RAD soliton or the depth of the RD soliton is obtained as A\ — r Sp J * • The velocities of both the 
RD and RAD solitons are proportional to the amplitude of the cw wave, that is, v\ = —2 ap. In particular, 
for oTm(si) > 0, we use = arccos [ ( 2 p^un(^|)[+i)^ ] t0 characterize the darkness of the RD soliton. 
If cci = | (i.e., |Im(si)| = ^), u\ represents a rational black soliton; while for 0 < cci < | (i.e., 
|Im(si)| / 0,7^;), it describes a rational gray soliton. 


7 












Figure 2: Degenerate two-soliton interactions via the solution in Eq. (fl4l) : (a) Segmented RAD soliton with p = 
<j) = 0, a = —1, si = 1. (b) Segmented RAD soliton with p = <fi = 0, a = —1, si = y/2i. 


Second, we derive the asymptotic expression of the solution in Eq. (fl4l) along the line x — 2 apz ^ 0 as 
z\ -A oo as follows: 


-Mif :=pe Tlp2z+ * 


M s 


i n 12 2 

Kl =p 


1 + 


Ap{x — 2 crpz) — 4 if* + 4 
2p[x -2 apz) - K* 

8 (crplm(si) + 1) 

\2p(x — 2 apz) — K| 2 


(17a) 

(17b) 


Here, the intensity \u\\ 2 can also display the RD and RAD soliton profiles which are associated with 1 + 
aplm(si) < 0 and l+crplm(si) > 0, respectively. In this case, both the RD and RAD solitons are localized 
along the line x — 2 apz + Re(si) = 0, and thus their velocities are given by V 2 = 2 ap. The height of \u\\ 2 
for the RAD soliton or the depth of \u\\ 2 for the RD soliton is equal to A 2 = ' Similarly, 

for 1 + crplm(si) < 0, we can use CZ 2 = arccos r 8 (pl Im ( 5 i)l~ 1 ) 1 to characterize the darkness of the RD 

L ( 2 p|Im(si)| — 1 ) J 

soliton. When 002 = f (i.e., |Im(si)| = 2v ^~ 1 ), is a rational black soliton; while for 0 < u>2 < f (i.e., 
|Im(si)| / 2 ^ 2 p l ? if represents a rational gray soliton. 

The above asymptotic analysis implies that the solution in Eq. (fl4l) can describe the elastic interactions 
of rational solitons in the sense that two interacting solitons retain their individual shapes, intensities, and 
velocities as z -A ± 00 . However, different from the standard elastic interaction in the NLS model, each 
interacting soliton experiences no phase shift upon an interaction. In general, the solution in Eq. (fl4l) exhibits 
three different types of elastic interactions between two rational solitons on a cw background, as shown in 
Fig. [I] The associated parametric conditions are given in the first three rows of Table Q] In particular, with 
Im(si) = 0, the asymptotic soliton u\ vanishes as z -A ± 00 , while u\ displays a RAD soliton profile [see 
Fig. |2(a)| . Similarly, for the degenerate case 1 + crplm(si) = 0, the only surviving asymptotic soliton is u\ 
and it takes the shape of the RAD type [see Fig. |2(b)] . In either of the two degenerate cases, one asymptotic 
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Table 1: Asymptotic patterns of the solution in Eq. (fl4l) under different parametric conditions. 


Parametric conditions 

Asymptotic soliton u\ 

Asymptotic soliton u\ 

crlm(si) < 0 , 1 + crplm(si) > 0 

RAD soliton 

RAD soliton 

(jlm(si) < 0, 1 + < 7 plm(si) < 0 

RAD soliton 

RD soliton 

(jlm(si) > 0, 1 + < 7 plm(si) > 0 

RD soliton 

RAD soliton 

Im(si) = 0 

Vanish 

RAD soliton 

1 + crplm(si) 0 

RAD soliton 

Vanish 


soliton disappears in the far-held region, but it still affects the other one in the near-held region, that is, the 
surviving soliton is segmented into two pieces at some hnite value of z. Therefore, such two degenerate 
cases of the solution in Eq. (fl4l) cannot be simply regarded as the conventional single soliton [36]. 

3.2 Second-order rational soliton solution 

With the truncation of Eqs. (l 8 al) and (l 8 bl) at j > 2, we can further obtain a series of higher-order rational 
soliton solutions. For j = 2, the second-order rational solution can be obtained as 


with 


u 2 = P e 2 ' lp2z+i *+ 2 


tr 

y(0,3) 

_(0,1) 

~9i 

rr 

ir 

fr 

rr 

-sr 

tr 

/r 

yr(0,3) 

r 3) 

-(0,1) 

-(i,i) 

~9 1 

~9 1 

=(0,3) 

=(1,3) 

~9 i 

~9i 


r ] 

(0,1) 

Pi 

y(2,3) 


-(2,1) 

“Pi 

/r 

-(2,3) 

“Pi 

fl° ,3) 


(o,i) ( 1 , 1 ) ’ 

Pi ~9i 

(0,3) (1,3) 

Pi ~9i 

j(0,l) _j( 1,1) 

J( 0,3) _jF( 1,3) 


(18) 


f^=2V2ipe^ C , fr=ia Pf r\ = 

% (cpC + i)> = i ppgf’^ 


(°4) = 2V2e- ip2z ~% _l_ 5 \ J 1 - 1 ) - 


9i 


e (p,3) _ ipe I/,2z+ 2 (4 p 2 ( 3 - 3C - 12 w) 

h ~ 


9? i) =-P 2 9 r\ 


3y/2 


,(i,3) _ P 2 cre ip2z+1 2 (4p 2 C 3 - 15C - 12 w) 
h ~ 3\/2 

(4p 2 C 3 -27C- 12 ro) 


,(2,3) _ ip e 

J l — 


3„i P 2 z+'4 


3\/2 


Pi 


(0,3) _ _£ 


ip2z ~'i (4(7/9 3 C 3 + 12i / o 2 C 2 - 15apC - 12 a pm - 3 i) 


3\/2 


(19a) 

(19b) 

(19c) 

(19d) 

(19e) 

(19f) 
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(1.3) a P e 

9i = - 

(2.3) P 2q 

9i =- 


-ip 2 z~'4 


• 2 

■'p z ~f 


(—4i ap 3 ( 3 + 12 /? 2 £ 2 + TJiapC, + 12i apw — 15) 
3\/2 

(4 ap 3 ^ 3 + 12 i/? 2 £ 2 — 39<jp( — 12 up w — 27 i) 

371 ’ 


(19g) 

(19h) 


where 


_ *(m—l,j) 

1 


(— x, z) (1 < m,j < 3), ( — x + 2 crpz + s i, 


/ \ —IT? 

w = 2 apz + s 2 , «i and 52 are two complex parameters. 

For the second-order rational solution in Eq. (fl8l) . one can find that it exhibits the finite-amplitude 
localized wave structures if and only if s± and 52 satisfy the nonsingular condition 


[2plm(si) + a] [2p 2 Im(si) 3 + 3palm(si) 2 — 3Im(s 2 )] < 0. (20) 

However, in contrast to the exponential soliton solutions of Eq. <[T|) |36], the second-order rational solution 
does not describe the interactions among larger numbers of RD and RAD solitons such as u\ or u\. Our 
asymptotic analysis reveals that there are only two asymptotic expressions when \x — 2 crpz\ —>> 00 and 
\x + 2 crpz\ —>>00, which are respectively given as 


u 2 -7 u\ := pe 2ip2z+ ^ 


u 2 u\ := e 2ip2z+ ^ 


_ 6 i cr p(€ + Sl ) 2 _- 

2 p 2 (£ + si) 3 + 3 i ap(£ + 5i) 2 + 3(2 apz + s 2 ) J ’ 

6 ia + 12/9(77 — s i ) — 6 icr/ 9 2 (77 — sl ) 2 
2 p 2 {rj — s^) 3 + 3iap(rj — s^) 2 — 3(2 apz + s^)- 5 


( 21 ) 

( 22 ) 


with £ and 77 defined below Eq. (fl4l) . By a numerical comparison, it can be verified that the solution in 
Eq. (fl8l) agrees very well with Eqs. (1211) and (l22l) at large values of z. Through the qualitative analysis of 
Eqs. (fl8l) . (l2TT) and (l22l) and the results in Figs. [3]-[71 we find that the second-order rational solution has the 
following properties: 


(i) In the near-field region \z\ <C 00, the intensity profiles of u\ and u\ display various combined-peak- 
valley soliton structures, which contain one or two peak(s) and valley(s) and vary with the evolution of 
z. As | z | —>► 00, either u\ or will slowly evolve into a RAD soliton or a RD soliton if Im(si) / 0 


and 1 + aplm(si) / 0 [see Figs. |3(b)H5(b)] l, and u\ 
aplm(si) = 0) [see Figs. |6(b)] and |7(b) [. 


(or U 2 ) will vanish if Im(si) = 0 (or 1 + 


(ii) The combined-peak-valley solitons u\ and interact elastically in the near-field regions [as shown 
in Figs. |3(a)l ~ [7(a)B , that is, the shapes, amplitudes, velocities, and phases of two interacting solitons 
remain the same after their mutual interactions. 


(iii) The parameter s\ determines whether u\ and u\ will eventually vanish or evolve into a RD or RAD 
soliton as | z | —» 00, while s 2 only affects the peak height and the valley depth for each interacting 
soliton in the near-field regions. 

(iv) According to the ultimate states of u\ and u\ as \z\ —>> 00, the second-order rational solution can 
exhibit three different types of elastic interactions (Figs.[3]-[5]) and two degenerate cases (Figs.[6]and[7]). 
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Moreover, the parametric conditions of the five cases are the same as those listed in Table \T\ for the 
first-order rational solution (fl4l) . 



(a) (b) 

Figure 3: First type of soliton interaction via the second-order rational solution in Eq. (fl8l ) with p = 0.5, = 0, 

a = — 1, si = 1 + 1.5 i, S 2 = 1 + i. (a) Rational soliton interaction in the near-field region, (b) Transverse plots of 
two asymptotic solitons at different values of z. 



Figure 4: Second type of soliton interaction via the second-order rational solution in Eq. (fl8l ) with p = 1, </> = 0, 
a = —1, si = 10 + 5 i, S 2 = —4 + 100 i. (a) Rational soliton interaction in the near-field region, (b) Transverse plots 
of two asymptotic solitons at different values of z. 
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(a) (b) 

Figure 5: Third type of soliton interaction via the second-order rational solution in Eq. (fl~8l) with p = 1, </> = 0, 
a = — 1, si = — i, S 2 = 1 — 12 i. (a) Rational soliton interaction in the near-field region, (b) Transverse plots of two 
asymptotic solitons at different values of z. 
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(a) (b) 

Figure 6: First degenerate case of soliton interaction via the second-order rational solution in Eq. (fl8l) with p — 1, 
(j) = 0, a = — 1, si = — 1, 52 = 1 — i. (a) Rational soliton interaction in the near-field region, (b) Transverse plots of 
two asymptotic solitons at different values of z. 
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Figure 7: Second degenerate case of soliton interaction via the second-order rational solution in Eq. (f]~ 8 l) with p = 0.5, 
(j) = 0, a = — 1, si = 1 + 2 i, 82 = 0.2 i. (a) Rational soliton interaction in the near-field region, (b) Transverse plots 
of two asymptotic solitons at different values of z. 


3.3 Stability analysis via numeric simulation 

On one hand, we use the time-splitting Fourier method to study the stability of the rational soliton 
solutions with respect to finite initial perturbations. We choose the solutions in Eqs. (fl4l) and (fl8l) at z = — 20 
as the initial values, and examine two types of initial perturbations: (i) small white noise on the initial values 
and (ii) small perturbation on the initial amplitude. In Figs. |8(a)| and |9(a)[ the numerical simulations present 
the stable evolution of the first- and second-order rational soliton solutions with the same parameters as 
those in Figs. |l(c)| and [5l respectively. When a white noise with a maximal value of 0.1 is added to the 
initial values, it can be seen from Figs. |8(b)1 and |9(b)l that the propagations of the first- and second-order 
rational soliton solutions are affected very little by the white noise. If the amplitudes of the initial values are 
amplified by 10%, Figs. |8(c)] and |9(c)| show that the magnitudes of two interacting solitons are enhanced but 
the soliton shapes are maintained very well. 

On the other hand, we note that the exact solution of Eq. (OD) can form a VT -symmetric self-induced 
potential that maintains the stable localized soliton structures. Here, we are concerned with whether the 
VT -symmetry breaking of the self-induced potential will lead to the instability of rational solitons. First, 
we consider that the initial values of u(x,z) and u*(—x,z) have the same shift in the ^-coordinate, i.e., 
u(x,z) —)• u(x — xo,z), u*(—x,z ) u*(—x + xo,z). In this case, there is no instability occurring in 

the evolution of the first- and second-order rational soliton solutions [see Figs. |8(d)1 and |9(d)] , although 
the symmetric center of the self-induced potential is shifted to xq- The reason lies in the fact that Eq. (OQ) 
remains invariant under the coordinate transformation x x + xo. Second, for the opposite shift in the 
x-axis u(x,z) —> u(x — xo, z), u*(—x,z ) u*(—x — xo, z), the self-induced potential cannot keep the 

PT-symmetry with respect to any point of x. Hence, an obvious instability will appear for both the first- 
and second-order rational soliton solutions, as displayed in Figs. |8(e)| and |9(e)| As the value of xq increases, 
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the soliton solutions become more unstable and the localized structures are finally destroyed. This is similar 
to the situation for the exponential soliton solutions of Eq. (00) [36]. 





(a) 


(b) 


(c) 



(d) (e) 

Figure 8: Numerical simulations of the first-order rational soliton solution, (a) The initial value is the exact solution 
in Eq. (fl4l) at 2 = —20 with p = 1, </> = 0, a = 1, si = 1 + i. (b) A white noise with a maximal value of 0.1 is added 
to the initial value, (c) The amplitude of the initial value is amplified by 10%. (d) u(x, —20) and u*(—x, —20) have 
the same shift xq = 6. (e) u(x, —20) and u*(—x, —20) have the opposite shift xq = 0.3. 
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(a) 


(b) 


(c) 



(d) (e) 

Figure 9: Numerical simulations of the second-order rational soliton solution, (a) The initial value is the exact solution 
in Eq. (fl8ll at z = —20 with p = 1, (j> = 0, a = —1, si = —i, S 2 = 1 — 12 i. (b) A white noise with a maximal 
value of 0.1 is added to the initial value, (c) The amplitude of the initial value is amplified by 10%. (d) u(x, —20) and 
u*(—x, —20) have the same shift xq = 9. (e) u(x. —20) and u*(—x, —20) have the opposite shift xq = 0.2. 
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4 Conclusions 


We have studied the rational nonlinear localized wave phenomena on the cw background for the VT - 
symmetric nonlocal NLS model with the defocusing-type nonlinearity. By the generalized DT, we have 
derived the nonsingular rational soliton solutions starting from a cw solution. For the first-order rational 
solution, we have revealed the elastic RAD-RAD, RD-RAD, and RAD-RD soliton interactions, in which 
there is no phase shift for interacting solitons. Meanwhile, we have discussed the degenerate case in which 
only one RD or RAD soliton survives. We have found that the second-order rational soliton solution does 
not exhibit the elastic interactions among larger numbers of fundamental rational solitons such as u\ in 
Eq. (116al) or u in Eq. (117al) . but the two interacting solitons display a rich variety of combined-peak-valley 
structures in the near-held regions, and each of them eventually vanishes or evolves into a RD or RAD 
soliton as \z\ oo. Also, we have numerically analyzed the stability of the first- and second-order rational 
soliton solutions. The results show that the soliton structures are stable with the addition of a white noise to 
the initial value and under small perturbations of the amplitude of the initial value, but the stability will be 
destroyed if the self-induced potential loses the VT symmetry with respect to any point of x. 
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